Measurements of identified particles at intermediate transverse momentum in the 
STAR experiment from Au+Au collisions at y^s/vjv = 200 GeV 
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Data for Au+Au collisions at y^sjvjv ~ 200 GeV are analyzed to determine the ratios of identified 
hadrons (tt, K, p, A) as functions of collision centrality and transverse momentum (pr). We find that 
ratios of anti-baryon to baryon yields are independent of pr up to 5 GeV /c, a result inconsistent with 
results of theoretical pQCD calculations that predict a decrease due to a stronger contribution from 
valence quark scattering. For both strange and non-strange species, strong baryon enhancements 
relative to meson yields are observed as a function of collision centrality in the intermediate pr 
region, leading to p/vr and A/K ratios greater than unity. The increased px range offered by 
the A/Kg ratio allows a test of the applicability of various models developed for the intermediate 
Pt region. The physics implications of these measurements are discussed with regard to difi'erent 
theoretical models. 

PACS numbers: 25.75.Dw, 25.75.-q 



I. INTRODUCTION 

The Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL) is dedicated to 
the production and study of strongly interacting matter 
under extreme conditions of temperature and density, 
formed when nuclei at very high energies collide. Lat- 
tice QCD predicts for such extreme conditions that the 
strongly interacting matter will undergo a phase transi- 
tion leading to a system where the degrees of freedom are 
partonic Experimentally, however, it is not pos- 

sible to observe such a partonic system, but rather one 
observes a system that has evolved from the partonic sys- 
tem and whose degrees of freedom are hadronic. These 
hadronic data must then be used to infer information 
about the original partonic system and its evolution. 

It has been shown that a hydrodynamic description of 
relativistic heavy-ion collisions reproduces the main char- 
acteristics of soft particle production in the region of rel- 
atively low transverse momentum (~ pT < 2 GeV /c) 
At high momentum transfer px > 5 GeV/c), parton- 
parton scattering and subsequent string fragmentation 
should dominate particle production, and one expects 
theoretical pQCD calculations to describe the data aris- 
ing from hard processes. The intermediate transverse 
momentum region (~ 2 < px < 5 GeV/c) sits tacitly 
between these two extremes, and here one might expect 
data to reflect an interplay between both mechanisms 
rather than a simple superposition; not only could par- 
ton fragmentation be influenced by the surrounding mat- 
ter at high Pt 0, but non-perturbative effects might 
be observed just above the region where the hydrody- 
namical description can successfully reproduce the mea- 
surements 0. Indeed, it has been suggested that at 
top RHIC energies, specific non-perturbative mechanisms 
might exist at intermediate pr Q, H, |^ ^3 ■ ^^i^ paper 
we present new data which extend to the upper bound of 
the intermediate px region and make comparisons of the 
data with model predictions which incorporate a variety 
of mechanisms. 

The paper is outlined as follows. Section II provides a 
brief description of the experiment which took data for 
AuH-Au collisions at ^/SNN = 200 GeV. The details of 
the analysis, using data from a ring imaging Cherenkov 
detector and a weak decay topological method, are 



presented in Section III. The analysis allows charged 
particle identification (p, p, K"*", K~, ir'^ and 7r~) 
and neutral hadron reconstruction (A and Kg) in 
the intermediate px region. The results (px spectra, 
anti-baryon over baryon ratios and baryon over meson 
ratios) are presented in Section IV. Comparisons of data 
with available model predictions are made in Section 
V. Finally, Section VI contains a summary and some 
conclusions. 



II. EXPERIMENT 

During Run 2 at RHIC (years 2001-02) data were 
taken using the STAR detector [ill for Au-hAu col- 
lisions at 200 GeV per nucleon pair. For tracking 
and momentum determination in the fiducial volume, 
star's magnetic field was 0.5 T. For charged particle 
identification (PID), both the full azimuthal coverage of 
the Time Projection Chamber (TPC) ^1 and a Ring 
Imaging Cherenkov detector (RICH) [l3 | were used. 
The latter has smaller acceptance than the TPC but 
has very good particle identification capabilities in the 
intermediate momentum range. 

Two triggers were used to record data. The first, the 
minimum bias trigger, required a coincidence among two 
Zero-Degree-Calorimeters (ZDCs) located ~ 18 m up- 
stream and ~ 18 m downstream from the center of STAR 
and very close to the beam axis (8 < 2 mrad) 0|. The 
ZDCs detect spectator neutrons — neutrons in the collid- 
ing Au nuclei that do not participate in the Au+Au colli- 
sions. The second trigger, a charged particle multiplicity 
or centrality trigger, required, in addition to ZDC coin- 
cidences, hits in STAR'S Central Trigger Barrel (CTB). 
The CTB is an array of scintillator slats that surround 
the TPC at a radius of 2 m. The CTB signal is propor- 
tional to the charged particle multiplicity in the pseudo- 
rapidity interval \ri\ < 1. Use of the CTB in the trigger 
allowed the selection of events according to their cen- 
trahty class, up to approximately the top [0-5]% of the 
inelastic cross-section In addition, for the highest 

luminosities and data taking rates, the average number 
of high transverse momentum particles within an event 
was enhanced through fast online reconstruction for par- 
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tides hitting the RICH detector [l^ . 



III. ANALYSIS 

Data for 1.6 x 10^ Au+Au coUisions were recorded us- 
ing the minimum bias trigger, while the centrality trigger 
was used to record another 1.6 x 10^ events. Tracks were 
reconstructed in the TPC after an off-line selection on 
the position of the primary vertex along the beam axis 
(|z| < 25 cm, where z = is the center of STAR) this was 
required to avoid selection bias when using ZDC timing 
only. This vertex cut also limited acceptance bias due to 
the CTB response for selecting high multiplicity events 
with the central trigger. The central and minimum bias 
samples were divided into centrality classes. Data for the 
four classes [0-5%], [20-40%], [40-60%] and [60-80%] are 
presented and discussed in this paper. The most cen- 
tral class contained 800 k events, while each of the other 
classes contained 325 k events. Use of the RICH detector 
provided a smaller subset of events, approximately 80% 
of each sample. 



A. Pion, kaon and proton identification with the 
RICH detector 

The STAR Ring Imaging Cherenkov detector (RICH) 
is 1 m^ in area and was positioned at central rapidity a 
radial distance of 2 m from the beamline. It's azimuthal 
coverage was Acj) < 20° which corresponds to approxi- 
mately 2% of the total TPC acceptance [iTl j. 
The momentum of each charged particle was determined 
from the curvature of its reconstructed helical track in the 
TPC. After reconstruction in the TPC fiducial volume, 
some tracks were extrapolated to the RICH. The corre- 
sponding minimum ionizing particle (mip) first crosses 
the Cherenkov photon radiating medium. The mip ra- 
diates photons while going through this 1 cm thick liq- 
uid (minimum for a normal incident angle). A 0.5 cm 
quartz window with very good transmission for created 
UV-photons [T^ separates the radiating medium from a 
proximity gap of 8 cm. Both mips and induced photons 
reach the photo-detector consisting of a multi-wire pro- 
portional chamber with a Csl layer deposited on the pad 
cathode surface for photo-conversion. Rings are recon- 
structed from the photo-electron cluster locations due 
to the Cherenkov light and the central cluster due to 
the MIP. For each particle, a Cherenkov angle is com- 
puted using the ring and track information at the radia- 
tor; then the relativistic velocity, /?, is obtained using the 
Cherenkov angle once the momentum is known. 

Primary tracks that are reconstructed in the TPC and 
within the acceptance of the RICH (jr^j < 0.2) are se- 
lected with a distance of closest approach (DCA) to the 
primary vertex < 0.75 cm. This selection increases the 
probability that tracks originate from the interaction ver- 
tex (primary tracks) rather than from secondary decay 



TPC Geometrical Selections 




Value 


Pt 


> 


1.5 GeV/c 


\v\ 


< 


0.20 


DCA to primary vertex 


< 


0.75 cm 


Fit hits 


> 


20 


Last hit radial distance ^ 


> 


140 cm 


Ring Imaging Selections 




Value 


MIP cluster charge 


> 


150 ADC 


Ring photons 


> 


1 


Incident angle 


< 


17° 


Residual 


< 


0.5 cm 



^ maximum number of hits is 45; 



^ outer radius of the TPC is 200 cm. 



TABLE I: Values of the geometrical cut variables used to 
select the TPC tracks entering the RICH detector and values 
of the cut parameters for the RICH analysis. 

products or particles produced via interactions in detec- 
tor material. Other selections are required in order to 
have sufficient track quality before extrapolation to the 
RICH detector: a minimum of 20 hits in the TPC (out of 
a maximum of 45) with the most external one at a radial 
distance of at least 140 cm. 




Photon Cluster Charge (ADC) MIP Cluster Charge (ADC) 

FIG. 1: Distributions of cluster charge for all the detected 
photons (left panel) and for tlie mip after pattern recognition 
(right panel). The final value for the mip selection is indicated 
by an arrow. 



The energy deposited on the pad plane, larger for a 
MIP, can be used to distinguish charged particles from 
photons. A cluster charge of more than 150 ADC counts, 
as shown in Fig.^ identifies the mip without ambiguity. 
Its centroid defines the center for the ring recognition 
algorithm. It is possible that a ring can be determined 
by this center and a single photon only (cluster charge 
< 150 ADCs). Ring identification is improved by re- 
quiring a small incident angle between the detector axis 
and the direction of the track (< 17°). A comparison 
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of the extrapolated TPC track and the center of the 
cluster from the mip leads to residual distributions which 
depend on the position of the interaction vertex along 
the beam axis. The related widths are 0.17 cm and 
0.20 cm in the local x and y directions of the pad plane 
respectively (i.e. parallel to the azimuthal direction and 
to the beam direction). These distributions arc shown 
in Fig. 121 The TPC and RICH selection criteria are 
summarized in Table 




-40 -20 20 40 -1 -0.5 0.5 1 

Primary Vertex z co-ord. (cm) Extrapolation Residuals (cm) 



FIG. 2: Left panel: The distribution of the primary vertex 
location along the beam axis before (solid) and after (dashed) 
centrality and trigger selections. Right panel: Residuals of 
TPC extrapolated tracks and center of mip cluster. Both 
figures are for the most central event sample. The small offset 
for the residuals along the x co-ordinate of the RICH plane is 
due to a systematic shift in the mean position of the primary 
vertex. 



The Cherenkov angle versus momentum for recon- 
structed charged particles is presented in Fig. 13 Lines 
correspond to the expected angles for pions, kaons and 
protons for an index of refraction of n = 1.29039 -By 
binning in transverse momentum, histograms such as the 
one shown in the insert of Fig. |2| enable particle identifi- 
cation with the RICH as discussed below. 



Identification of mesons is possible up to pt = 3 GeV/c 
and up to pt = 4.5 GeV/c for protons and anti-protons 
once the cuts in Tableware applied. The momentum of 
each particle at the primary vertex is different from the 
momentum at the RICH radiator. In order to perform 
a three-Gaussian fit to the Cherenkov angle distribution 
in each pT interval and extract the raw yields, the mo- 
mentum at the radiator level is needed. Therefore, the 
mean energy loss in the TPC and its outer field cage for 
the three different species (pion, kaon, and proton), as a 
function of momentum, is determined via Monte Carlo 
simulation. After full propagation with GEANT and 




Momentum (GeV/c) 



FIG. 3: (Color online) Cherenkov angle vs. momentum (mo- 
mentum obtained with the TPC). The inset shows the pro- 
jection of the Cherenkov angle for a small pr interval. 

reconstruction, Gaussian fit parameters for the Monte 
Carlo spectra are compared to those extracted from real 
data. 

We find very good agreement at low transverse momenta 
{pt < 2 GeV/c), i.e. where the three peaks are suffi- 
ciently separated for unambiguous determination of the 
mean and width of each Gaussian fit to the real data. 
With increasing px, the Gaussians start merging, and 
several fit strategies arc used in order to extract the raw 
yields: i) require that the three-Gaussian integral equals 
the total number of projected tracks; ii) use the means 
and/or the widths extracted from simulation to add more 
contraints; iii) when kaons and pions arc difficult to dis- 
entangle (for Pt > 2.6 GeV/c), fix the kaon yields with 
the measured Kg yields, presented later in this paper, as 
a unique constraint for obtaining the raw yields of pro- 
tons, anti-protons and pions. An example of a fit where 
the means of the Gaussians are fixed using simulation 
parameters is shown in Fig. 0] The arrows indicate the 
positions of the Gaussian means of the fits to the simula- 
tion data. The errors shown in Fig.^arc quadratic sums 
of statistical errors and systematic uncertainties due to 
the fitting. The overall systematic uncertainties for the 
extracted raw yields are obtained by changing the fitting 
techniques, varying the number of free parameters (leav- 
ing free and fixing means, widths, or requiring that the 
overall integral is equal to the number of counts) , as well 
as varying the cuts of Table n 



B. Topological identification of A and with the 
TPC. 

The strange hadrons, A and Kg (so called "VOs"), are 
reconstructed in the Time Projection Chamber at mid- 
rapidity (|y| < 1) via measurements of their charged de- 
cay products (daughters) from weak decays: A(A) — > 



6 



100 


2.6 < p-p (GeV/c) < 2.8 


A 






- □ Proton = 663±32 1 1\ 


1 






Z 1 Kann - ^70+?? / i 






50 


- Pion = 550±27 1 ^1 








Global fit / 1 1 








\ ,,,,,, , X 






3 0.4 0.5 0.6 


0.7 



Cherenkov Angle (rad) 

FIG. 4: (Color online) Extracted Cherenkov angle distribu- 
tion for particles with a positive charge and 2.6 < pr < 
2.8 GeV/c. This pr range is close to the limit where mesons 
cannot be identified using the integrated fit technique. The 
data used are those for the most central collisions. Using a 
constraint on the total number of counts gives a x'^/^DF — 
3.8. 



p{p) + 7r-(7r+) (64%) and K*^ ^ Tr+ + tt" (69%). The 
strange hadrons are reconstructed by pairing tracks of 
opposite sign and a small distance of closest approach 
(DCA) , while ensuring that the momentum vector of the 
assumed neutral parent is on a line through the interac- 
tion vertex. However, due to the large particle multiplic- 
ities in relativistic heavy-ion collisions, particle spectra 
produced in this way suffer from very large combinato- 
rial backgrounds, especially at low pT, owing to random 
crossings of uncorrelated tracks predominantly close to 
the interaction vertex. In order to minimize this back- 
ground, geometrical selections (i.e., cuts) are applied to 
the data which reduce the background while keeping as 
much of the signal as possible. These selections include 
imposing a minimum decay length for the VO, which re- 
moves all crossings of primary tracks at the interaction 
vertex which would appear as VOs with zero decay length, 
as well as cuts on the minimum DCA of the daughter 
particles to the interaction vertex. These selections are 
effective as the momentum vectors from true daughters 
do not lie along a line through the interaction vertex 
owing to the momentum kick given to these particles in 
the decay (larger for A'^ compared to A and A due to 
the differences in mass between the parent and sum of 
the daughter particles). However, this DCA sensitivity 
is found to be strongly dependent upon the pr of the 
parent particle, hence different selections were utilized as 
a function of the parent px- The values used for these 
selections for the three different particles are given in Ta- 
blem The raw mass spectra for the different hadrons arc 
shown in Fig. |S1 for a typical low transverse momentum 
bin where the backgrounds are the most significant. 
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FIG. 5: The relative raw mass spectra for the Kg, A and A 
candidates, respectively, within a given low transverse mo- 
mentum bin. 

In order to extract raw yields, two different methods are 
employed. The first uses a fit function to describe both 
the signal (single-Gaussian or double-Gaussian) and the 
background on either side of the peak (third-order poly- 
nomial). For the example shown in Fig. O the single 
Gaussian method fails to fit the signal, both in the tail re- 
gions of the distribution and the magnitude of the peak; 
however, for the example shown, the double Gaussian 
method reproduces the data well with a small for the 
fit and does so in the majority of the other cases. Oc- 
casionally, however, this method fails, and the fit does 
not converge. To avoid such cases, another method of 
extracting the raw yields is used. This method involves 
the summing of the entries of relevant bins in the his- 
togram as shown in Fig. |S1 and is referred to as 'bin- 
counting' Three regions of the histogram are iden- 
tified: sidebands on each side of the mass peak (each of 
width 15 MeV taken to represent the sum of the back- 
ground under the peak) and the signal area itself (in this 
case, ±15 MeV of the mass of the Kg) which contains 
both signal and background. The raw yield is the differ- 
ence between the signal area and the background area. 
This method works best if the shape of the background 
is linear over the full range. If the background is small in 
comparison to the signal, then any small deviations from 
linearity of the background have negligible effect on the 
extracted yield. 



C. Correction Factors 

After the raw yields are extracted as a function of pr, 
they must be corrected for analysis inefficiencies in order 
to obtain final spectra. These inefficiencies are calculated 
by embedding Monte Carlo (MC) tracks into real data 
and determining the probability of finding these particles 
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Geometrical Selections A A 

DCA VO to primary vertex < 0.5 cm < 0.5 cm < 0.75 cm 

DCA neg daughter to primary vertex" > 2.5 cm : > 4.0 cm > 1.5 cm : > 0.0 cm > 3.5 cm : > 2.0 cm : > 0.6 cm 

DCA pos daughter to primary vertex" > 1.5 cm : > 0.0 cm > 2.5 cm : > 4.0 cm > 3.5 cm : > 2.0 cm : > 0.6 cm 

DCA daughters at secondary vertex < 0.8 cm < 0.8 cm < 0.8 cm 

Decay length > 6 cm > 6 cm > 5 cm 

"This cut is strongly dependent upon the transverse momentum (pt). For the A and A, two different pr ranges were used, 
Pt < 3 GeV/c, pt > 3 GeV/c. For the Kg, three pr ranges were used, pr < 1.0 GeV/c, 1.0 < pr < 3.5 GeV/c, pr > 3.5 

GeV/c. 

TABLE II: The values of the geometrical selection variables used to reconstruct the strange hadrons. 
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FIG. 6: The raw mass spectrum for the Kg for a given trans- 
verse momentum bin, with three different signal extraction 
methods indicated. 



number of reconstructed MC particles to the number of 
generated particles. As the majority of the background 
is dependent upon the track multiplicity in the event, it 
is found that the efficiency is also strongly dependent 
upon collision centrality. Care is taken to ensure that 
when the mixing is done, the multiplicity of the mixed 
event is not increased by more than 5% of its original 
value. This value is chosen so that the efficiency is not 
compromised while maximizing the available statistics. 
In addition to the TPC MC simulation, a RICH MC 
simulation is done. Here, the ADC values corresponding 
to both the mip and the photons generated as the 
particles pass through the radiator are generated and 
used in the efhciency calculation. This is done during 
the same reconstruction pass as for the TPC embedding. 
Examples of the centrality and px dependence of the 
efficiencies are shown in Figs. [7| and |S1 which show 
efficiencies for the most central and most peripheral cen- 
trality classes used in the analysis presented in this paper. 



through data analysis. This process is performed in 
STAR utilizing the GEANT code. The MC particles are 
generated with a user-defined realistic distribution in 
transverse momentum and rapidity space. As the real 
data have a finite interaction vertex distribution, care 
is taken to ensure that the generated particles have the 
same interaction vertex distribution as the real data. 
GEANT propagates the particles through the material 
present in STAR, and the resultant ionization deposited 
in the TPC is then effectively drifted to the read-out 
planes using a detailed simulation response program 
which provides ADC values from the cathode pads that 
are in the same format as that of the real data. The 
data from the simulation are then mixed with real data, 
and star's full tracking and reconstruction software 
is used on the mixed data to determine the efficiency 
for finding the MC particles. As the MC information 
is recorded separately for the embedded particles, it is 
easy to determine the efficiency of the STAR analysis 
software for finding the embedded particles. The same 
geometrical and multiplicity cuts are used for the 
mixed data as are used for the real data. The total 
efficiency, for each px bin, is simply the ratio of the 




Pj (GeV/c) 

FIG. 7: (Color online) Efficiencies as a function of pt and 
centrality for the Kg and A. The three ranges in the efficiency 
are due to changing of geometrical reconstruction cuts given 
in Table ^ chosen for optimizing the signal-to-noise ratios 
and maximizing the statistics at high pr- 
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FIG. 8: (Color online) Charged particle reconstruction and 
identification efficiencies as a function of pr (left panel) . Col- 
lision centrality effects, mainly due to TPC track reconstruc- 
tion, are shown for the anti-proton and the negative pion for 
the extreme centrality bins used (right panel). 

At relatively low transverse momenta a 10% difference is 
observed among charged mesons because kaons can de- 
cay before reaching the RICH fiducial volume (left panel 
of Fig.|Sl). Within statistical errors no significant differ- 
ence is observed between opposite charges of the same 
species, except for baryons, due to anti-proton absorp- 
tion in the detector material. Although not shown for 
clarity in the right panel of Fig. |H1 efficiencies for the 
intermediate centrahty bins ([20-40%] and [40-60%]) sit 
between the extremes within statistical errors. This de- 
pendence on centrality is not a feature of the RICH due 
to low occupancy, but rather comes from the track find- 
ing inefficiency in the TPC, as a track has to be found in 
the TPC before it can be extrapolated to the RICH. The 
number of photo-electrons in a cluster for an identified 
Cherenkov ring is a function of the relativistic velocity 
and therefore of the momentum of the particle at the ra- 
diator level; it increases with f3 and reaches up to 17 ± 2 
when P ~ 1 [Ullllllllll. The left panel of Fig.Hshows 
distributions of the number of photons for the limiting 
cases of protons which have a significantly lower velocity 
for a fixed momemtum bin. The lower limit for the mean 
number of photons (for 1.5 < pr < 1.6 (GeV/c)) is ap- 
proximately 5 for this species and is correctly reproduced 
by the simulation (dashed line). 



Due to energy loss along the path from the interaction 
vertex to the RICH pad plane, significant differences 
can be seen in the momentum of the same particle 
at different points along its path. The right panel of 
Fig. 13 presents distributions of the differences of these 
momenta for simulated protons at low momentum where 
the differences can be relatively significant: no difference 
is seen between the simulated momentum (MC) and the 
one reconstructed at the primary vertex (RC) within 
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FIG. 9: (left panel) A comparison of the number of photons 
for different pt bins for protons, (right panel) The energy 
loss effects on momentum determination: simulated (MC), 
reconstructed at the primary vertex (RC), reconstructed at 
the outermost point of the TPC (OUT), as well as at the 
RICH pad level (PAD). 

momentum determination uncertainties. Though a small 
shift is distinguishable for the momentum at the end of 
the TPC (OUT), the most significant one is at the RICH 
pad level (PAD) once the track has travelled through the 
TPC, the outer field cage and the Cherenkov radiator. 
The mean momentum shift resulting from the overall 
energy loss is a function of the track momentum. It 
is obtained from simulation and taken into account to 
determine the corrected mean Cherenkov angle for each 
Pt bin. 



IV. RESULTS 

Once all efficiencies and correction factors are calcu- 
lated (both as a function of pr and centrality), they are 
applied to the raw data, and final spectra are obtained 
as illustrated in Figs. ITUl [TTI andini 

Figure lTUl shows a comparison between the non-strange 
identified spectra from the RICH analysis and those de- 
termined from energy loss {dE/dx) in the gas volume of 
the TPC for the most central bin The extra reach 
of identified particles obtained using the RICH detector 
is significant. Figure lTTl shows the p and 7r~ spectra from 
the RICH for three of the centrality classes studied, while 
Figure [T^ shows the strange particle spectra for all four 
centralities used. 



The maximum extent in px for the A and Kg is not re- 
stricted by the identification method, but rather is lim- 
ited by statistics only. Therefore, the reach in px for the 
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FIG. 10: (Color online) Identified charged particle spectra 
from analysis of the RICH data (right of vertical line) ; to the 
left of the vertical line are identified charged particle spectra 
determined using energy loss (dE/dx) in the TPC 12^1 . 
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FIG. 11: (Color online) The p (full symbols) and tt (open 
symbols) spectra as a function of both centrality and pr- 
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FIG. 12: (Color online) The A (full symbols) and (open 
symbols) spectra as a function of both centrality and pr- 



[0-5%] centrality bin extends further than for the other 
centrahty bins as it includes data taken with the most 
central trigger, mixed with the most central portion of 
the data talcen with the minimum bias trigger. 

The results shown in Fig. ^] take into account both 
statistical and systematic errors (added in quadrature); 
their values arc smaller than the symbol sizes. For both 
particles, the systematic errors are dominated by sig- 
nal extraction uncertainties, which vary from 5(20)% at 
very low pr for K°{A), to 10(10)% at high pr- The 
other sources of systematic error include effects due to 
magnetic field polarity, effects due to differences between 
TPC East data and TPC West data, as well as effects due 
to the choice of geometrical cuts. We estimate that these 
effects contribute an additional 5% to the uncertainty for 
both the and A. 

Neither the anti-proton nor A hyperon spectra pre- 
sented in this paper are corrected for feed-down from 
weak decays. The dominant channels were studied for 
each case and found to be 15-20% for the anti-protons, 
10% for the A, and independent of transverse momentum. 

The p/p and A/A ratios are plotted in Fig. E| as a 
function of pr for the [0-5%] centrality bin. Both sta- 
tistical and systematic errors are shown. For both cases, 
raw particle yields are used because of the charge sym- 
metric nature of both the applied cuts for particle re- 
construction and the tracks. The ratios are corrected for 
absorption in the detector material using a GEANT sim- 
ulation; the largest absorption correction occurs for the 
anti-baryon. It is expected that at large px, where par- 
ticle production is dominated by jet fragmentation, the 
anti-baryon/baryon ratio (B/B) will start to decrease. 
This is because the leading baryon from a quark jet is 
more likely to be a particle rather than an anti-particle, 
whereas this asymmetry does not occur in a gluon jet. 
For the most central collisions, this trend is not observed 
in either the p/p ratio out to pr = 4.5 GeV/c or the 
A/A ratio, out to 5.5 GeV/c. The ratios for the other 
centrality bins are not plotted for clarity. However, as 
a function of centrality, neither ratio shows a decrease 
though the px coverage is less, due to lack of statistics. 
The integrated value of the ratio is constant for all cen- 
tralities. 



The baryon to meson ratios (B/M) can also be stud- 
ied with this data-set to investigate further the relative 
baryon enhancement first reported in the 130 GeV data. 
That work showed that the A/h~ ratio approached unity 
in the intermediate px region [l^ . as did the p/ir ra- 
tio [23]. This is contrary to what is observed for ele- 
mentary particle collisions where the ratio does not ex- 
ceed 0.2 in the corresponding momentum range for both 
quark and gluon jets This is as expected for jet 

fragmentation. 

FigurelTHshows the p/tt~ ratios, again for the 5% most 
central collisions, as well as for both the [20-40%] and 
[40-60%] centrality classes. The px range is limited to 
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FIG. 13: (Color online) The p/p (circles) and A/A (triangles) 
ratios as a function of pr- The statistical errors are indicated 
by the error bars, the systematic errors are indicated by the 
hash marks. 



3 GcV/c due to the need for unambiguous identification 
of the 7r~ and by statistics for the more peripheral sam- 
ples. The ratios clearly increase with increasing px and 
also exhibit centrality dependence; for a given px, the 
magnitude is larger for the more central collisions. How- 
ever, for the pt range studied here, there is no evidence 
for the ratio turning over or decreasing to the small values 
measured in jet production [2^. It should occur, how- 
ever, for high enough momenta when fragmentation is the 
dominant process. We note here that this ratio is con- 
sistent with the ratio presented by the PHENIX collabo- 
ration once the differences in centrality and treatment of 
feed down from weak decays are taken into account [2^ . 



possible to measure an approximation of the B/M ratio 
to higher values of by plotting the p/{h~ —p) ratio as a 
function of px, where the value h~ represents the sum of 
non-identified negatively charged hadrons {■k~^K~ and 
p) . As the identification of the particles is not required for 
the h~ data, it is possible to use the TPC data which has 
a larger acceptance than the RICH detector and hence 
the reach in px for the p/{h~ —p) ratio is dominated by 
the reach in px of the identified p. The p spectrum has 
been subtracted in the denominator because it is part 
of the measured h~ spectrum. This ratio is plotted in 
Fig. ^1 which shows that the reach in px extends to 
4.5 GeV/c. With the addition of the kaons in the denom- 
inator, the ratio shows an initial shape which is different 
from that of the p/tt" ratio, but the region of importance 
is for Px > i GeV/c where the trend of the data is for a 
slight decrease in the ratio. 
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FIG. 14: The p/tt ratio as a function of centrality and px 
for Au-|-Au collisions at ^/SNN ~ 200 GeV. 



As the identification of p and p using the STAR RICH 
detector extends well beyond that of the tt^, it is still 



FIG. 15: The p/{h —p) ratio as a function of px for Au+Au 
collisions at ^/snn ~ 200 GeV. 



In order to extend the reach in px even further, we 
use the strange particles which can be identified out to 
6 GeV/c in the most central bin (as illustrated in Fig. ll2|l . 
The A/ Kg ratio, a measure of the baryon/meson ratio in 
the strangeness sector, is plotted in Fig. ^| as a function 
of Px for the four centrality classes studied. The quoted 
error bars are statistical only as the calculated systematic 
errors for the A and Kg are correlated to an unknown 
degree, resulting in smaller systematic errors. With the 
increased px range for the B/M ratio, a definite turnover 
for strange particles is observed at pt ~ 3 GeV/c for all 
centrality classes. For px > -i GeV/c for each centrality 
bin the ratio decreases. At low transverse momenta the 
ratios appear to be independent of centrality (up to px ~ 
1.5 GeV/c). Beyond px ~ 1.5 GeV/c, at each px, the 
ratios increase with increasing centrality. 
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FIG. 16: (Color online) The A/Kg ratio as a function of 
centrality and pr for Au+Au collisions at ^/snn = 200 GeV. 



The A/ Kg ratio is greater than unity for the [0-5%] 
and [20-40%] data in the px range 1.5 < px < 5 GeV/c; 
in all cases it is larger than the B/M ratio measured for 
elementary collisions psf . 

For the first time, it is shown that the A/ Kg ratio 
(and hence the B/M ratio) does not keep rising with 
increasing px, but rather decreases for the region ^ 3 < 
Pt <5 GeV/c. 

The data suggest that at px 5 GeV/c the ratios for 
each centrality bin are roughly equal. This observation 
then suggests that the region ^ 2 < px < 5 GeV/c is the 
intermediate px region to use to test the ideas discussed 
in the introduction. 

In the following section, we compare a number of theo- 
retical models to the data and discuss their applicability. 
Each model uses a different method to describe particle 
production in the intermediate px region. 



tide production is not dominated by jets, but rather is 
governed by soft particle production and hydrodynam- 
ics 0,0. The 'Soft-|-Quench' model, however, is in good 
agreement with the data and predicts p/p and A/ A ratios 
which are flat up to about px ^ 5 GeV/c. This model 
links soft particle production at low px (where particle 
production is assumed to be proportional to e"^^/"^, T 
being the inverse slope parameter) with the higher px 
region dominated by a leading order pQCD calculation 
incorporating gluonic baryon junctions |28l |. The model 
also includes initial state multiple scattering (related to 
the "Cronin effect" 13), nuclear shadowing and par- 
tonic energy loss due to gluon bremsstrahlung as calcu- 
lated in the Gyulassy-Levai-Vitev (GLV) formalism [30l| . 
The third model prediction in Fig. 1171 is from a recombi- 
nation model [DU KE) which reproduces the p/p ratio 
well in the range where there is data The recombi- 
nation model itself is again governed by soft processes 
at low Px and at high px by jet fragmentation in ad- 
dition to jet quenching mechanisms; the predicted ratio 
starts to decrease for > 8 GeV/c. The model as- 
sumes that recombination processes (where two quarks 
combine into mesons and three combine to form baryons) 
are dominant at low px where the partonic spectrum is 
exponential in shape, while fragmentation becomes dom- 
inant in the high px region when the parton spectrum 
exhibits a power-law behavior. In this implementation, 
jet quenching mechanisms are also applied to the high 
Px region. Therefore, in this model, the intermediate px 
region (2 < pT < 5 GeV/c) is dominated by recombina- 
tion. 



Turning to the B/M ratios. Fig. El shows compar- 
isons between the p/'k^ ratio as a function of centrality 
and predictions from a variety of models. Besides re- 
sults from the models described above, predictions are 
also shown for a hydrodynamical model with a soft equa- 
tion of state l3ll , and from two additional recombination 
models HITc 



V. DISCUSSION 



The p/p and A/A ratios, presented earlier, are plot- 
ted in Fig. El a-s functions of px for the [0-5%] cen- 
trality bin, together with predictions of various the- 
oretical models. The main feature of both ratios is 
that they are essentially constant as a function of px- 
This observation contradicts predictions from pQCD in- 
spired models which show a decrease in the ratio over 
all vx 113, 113 ■ The prediction from one of these models 
([2g) is shown in Fig. El though the prediction is for 
Au-|-Au at y'sNN = 130 GeV. The reason for the de- 
crease is that in a quark jet, the leading baryon is more 
likely to be a particle than an anti-particle, whereas in 
a gluon jet, there is no such asymmetry. Although the 
pQCD prediction extends to below px = I GcV/c, it is 
not expected to be realistic in the low px region as par- 



Although hydrod yna mical models have been shown to 
work well at low px |31| , at the higher px (see left panel of 
Fig. I18|l . they cannot reproduce the data. In this model, 
a is proportional to an initial transverse seed velocity. 
The major discrepancy can be attributed to the pre- 
dicted 7r~ spectrum which deviates from the measured 
data at Pt ~ 1 GeV/c. Conversely, the 'Soft+Quench' 
model provides a qualitatively good description of the 
data, matching the observed rise in the data as well as 
the larger magnitude of the ratio for the more central 
data-set H^- For this model, there are predictions for 
two different geometrical scalings. Np^rt scaling refers 
to scaling with the number of participants, while Nf,m 
scaling is with respect to the number of binary collisions. 
The largest difference between these scalings occurs in 
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FIG. 17: (Color online) The p/p (upper plot) and A/A (lower 
plot) ratios as a function of pT, together with various model 
predictions (see text for full description). 



the more peripheral bin where the Npa^t scaling gives 
better agreement with the data. 

The 'Soft+Quench' model predicts the ratio's behav- 
ior to be due mainly to two effects. Firstly, pion pro- 
duction in central collisions is found experimentally to 
be lower than that in p + p collisions when scaled by 
the relevant number of binary collisions j26i | . This is at- 
tributed to quenching of jets due to high pT particles in- 
teracting with a dense gluonic medium. Secondly, baryon 
production is dominated at low to intermediate px by 
non-perturbative baryon junctions. The junctions them- 
selves are not found to be affected by the jet quenching, 
and hence the baryons scale with binary collisions from 
p+p data as expected. This scaling leads to an enhanced 
baryon/meson ratio at intermediate pt- It only decreases 
at higher momenta when perturbative processes become 
relevant for baryon production. The dependence of the 
ratio on ccntrality is naturally explained by the lack of 
jet quenching in the more peripheral collisions owing to 
the smaller, less dense medium that is created. 

On the right panel of Fig. ^1 three more model pre- 
dictions are compared to the experimental results. These 
models are similar in that they all subscribe to a recombi- 
nation (or coalescence) mechanism of dressed (i.e. mas- 
sive) quarks, but they have different implementations. 
In the following discussion, all references to partons per- 



taining to the coalescence models will refer to massive 
quarks. 

In each case, the models attribute particle produc- 
tion to the recombination mechanism when the par- 
ton transverse momentum spectra are exponential. At 
^ Pt > 5 GeV/c, where it is believed that the parton 
distribution is governed by a power-law, fragmentation 
dominates. Particle production at low pT (< 1 GcV/c) 
is governed by hydrodynamics. 

There are two methods used to evaluate the recombina- 
tion integral. The first is an analytical solution, though 
assumptions are made which cause a break down at low 
Pt- The second uses a Monte Carlo technique to evaluate 
the integral, negating the need for any approximations. 
Different approaches for the two methods are taken as 
described below. 

The first recombination model {DUKE - discussed 
earlier in the p/p section) uses the approximation tech- 
nique to do the recombination integral and allows for the 
recombination of uncorrelated thermal, valence quarks 
(i.e., quarks with an effective mass) only. The prediction 
is shown by the dashed line (right panel. Fig. I18|l . For 
the hard scattering part of the spectrum, the partons are 
taken from a mini-jet calculation which includes energy 
loss 0113. 

The second recombination model result shown in the 
right panel of Fig. [H by the sohd fine {TEXAS 0) 
uses the Monte Carlo method to do the integral. As 
well as allowing the recombination of thermal quarks, it 
also allows the recombination of thermal quarks with co- 
moving quarks coming from mini-jets 1^. This has the 
effect of increasing the ratio in the intermediate pt range. 

The third recombination model {OREGON [13) again 
allows for recombination of thermal quarks but differs 
again in how the higher pt regions are treated. Frag- 
mentation functions are not used in the intermediate pt 
region, but instead, mini-jets are allowed to form show- 
ers which then recombinc; this is the original use of re- 
combination models js^. The shower partons arc then 
allowed to recombinc with the thermal quarks. Results 
from two versions of this model are shown in Fig. 1181 
the dashed-dotted line represents the original prediction 
where the mass of the proton is neglected, and hence 
below Pt ^ 2 GeV/c, its assumptions become tenu- 
ous. The small dotted line represents the result of a later 
calculation where transverse mass is used instead of the 
proton mass [s^ ■ The original calculation reproduces the 
proton spectra well, whereas the later one is in agreement 
with the p/tt~ data. It is evident from Fig. El that the 
recombination models, though all showing a p/tt^ ratio 
which increases with pT up to about 3 GeV/c, have dif- 
ferent behaviors at lowp^. However, in common with the 
'Soft+Quench' model, they show a turnover to a falling 
ratio at Pt ~ 3 GeV/c, exactly where the data stop. 
Therefore, we use the ratio of strange baryons (A) to 
strange mesons {Kg) to investigate this turnover region. 
This comparison is made in Fig. 1191 
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FIG. 18: (Color online) The p/tt" ratio for the [0-5%], [20-40%] and [40-60%] most central collisions together with different 
model predictions for the most central bin (together with a prediction for a peripheral bin for the 'Soft+Quench' model). For 
the case of the DUKE model, the prediction is for p/tv". 
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FIG. 19: (Color online) The A/ Kg ratio as a function of centrality and pr for Au+Au collisions at y^sjvjv = 200 GeV, together 
with various model predictions. 



The left panel of Fig. ^| again shows the 
'Soft+Quench' model predictions, while the right panel 
only shows two predictions from the recombination mod- 
els {DUKE and TEXAS) as A/K^ predictions only exist 
for these models. OREGON model predictions for this 
ratio do not exist due to the lack of knowledge of the 
fragmentation function for hyperons at high px- Ad- 
ditionally, a prediction from a modified version of the 
TEXAS model is shown [s^. Note that the data have 
a slightly different centrality range than the range used 
by the models. Also, for the case of the 'Soft+Quench' 
model, the predictions are for the A/K^ ratio; this dif- 
ference should have minimal effect on the comparison be- 
tween data and prediction as the Kg yield is expected to 
be the mean of the K~^ and yields. It has been shown 
previously that the K~^ / K~ ratio is dependent upon the 
baryon chemical potential , and for Au+Au collisions 
at y/SNN = 200 GeV, this ratio differs from unity by 10% 
and is independent of pt |26j . 



The 'Soft+Quench' results are again in good agree- 
ment with the data for px > 5 GeV/c for the two cen- 
tralities shown. Below 3 GeV/c the model under-predicts 
the data, but reproduces the general shapes of the distri- 
butions. As for the recombination results (right panel of 
Fig. ll9|l . the DUKE model, which invokes the coalescence 
of only thermal partons, provides reasonable agreement 
with the data. The second set of model predictions are 
taken from the TEXAS model and show the result for the 
case where soft and hard partons are allowed to coalesce 
and the result from a modified version of the model in 
which only soft, thermal partons are allowed to coalesce. 
As expected, the difference between the two implementa- 
tions of coalescence manifests itself for ~ > 4 GeV/c, 
where the extra process leads to an enhancement above 
the data. In this representation, the versions of the mod- 
els which allow for the coalescence of thermal partons 
only agree better with the data. Due to the uncertainties 
in the fragmentation function for hyperons, the contribu- 
tion to the numerator for the A/ Kg ratios from fragmen- 



14 



tation is not included in this model, though this is only 
thought to be significant at higher px- The consequences 
of this is that at large pt, the ratio goes to zero instead 
of saturating at a value of approximately 0.1. 

Of the two classes of models discussed, the 
'Soft+Quench' model gives the best description of the 
data for both light and the strange hadrons for all ccn- 
tralities. However, in order to reproduce the data, novel 
baryon transport techniques must be added to the model. 
On the other hand, the recombination models provide an 
intuitive picture of particle production mechanisms for 
the intermediate pt range. They also reproduce other 
aspects of the data, such as the scaling of the elliptic 
flow with valence quark number, as measured by the 
STAR experiment |2a|- Though the different implemen- 
tations describe the central data, it would be instructive 
to see results of calculations for more peripheral collisions 
taking into account the full density profile of the cre- 
ated medium (39l] . In spite of the agreement with data, 
the current implementations of the recombination mech- 
anism can be improved. One of the main concerns is that 
entropy is not conserved due to the number of 2 ^ 1 and 
3 — !■ 1 processes which occur. This particular issue could 
be addressed by taking resonances into account, such as 
the p meson which decays into two pions. 

VI. SUMMARY AND CONCLUSIONS 

We have analysed -y/ijvw = 200 GeV Au+Au data 
taken with the STAR detector and its TPC and Ring 
Imaging Chcrenkov Detector (RICH) and extracted iden- 
tified particle spectra as a function of centrality out to 
high transverse momentum (px ~ 6 GeV/c). The p/p 
and A/ A ratios are both constant within errors as a func- 
tion of in the measured range (up to 5 GeV/c), which 
is inconsistent with results from pQCD calculations up 
to 4 GeV/c. At higher px the interpretation is incon- 
clusive owing to the larger error bars on the data. We 
have shown also that the p/tt" ratio rises with increas- 
ing pt up to ~ 3 GeV/c, while the A/ Kg ratio increase 



up to ^ 3 GeV/c and then turns over and starts to de- 
crease. In both cases, the ratio exceeds unity, meaning 
that there is a baryon dominance over meson yields in 
the intermediate pT region (~ 2-5 GeV/c). This is the 
first measurement to show a decrease at higher px and 
allows us to place an upper bound on intermediate pt 
region of approximately 6 GeV /c. 

We also compared the data with predictions from a 
variety of theoretical models. Even though the recombi- 
nation/coalescence models represent an intuitive physical 
process and can reproduce well the elliptic flow [4d| . 
when more comprehensive comparisons are made to 
ratios of particle spectra, significant differences are seen. 
Conversely, it was found that the 'Soft-|-Quench' model 
is better at reproducing both the pr and centrality 
dependence of the ratios, though this model contains 
more exotic mechanisms of baryon production (e.g. 
baryon junctions). 
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